Counting statistics for arbitrary cycles in quantum pumps 
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Statistics of charge transport in an adiabatic pump are determined by the dynamics of the scattering 
matrix S(t). We show that, up to an integer offset, the statistics depend only on the corresponding 
path N(i) = S^agS in the coset space (the sphere for a single channel). For a general loop S(t) we 
solve for the noise-minimizing pumping strategy. The average current is given by the area enclosed 
by N(t) in the coset space; its minimal noise by the area of a minimal surface (soap film) spanned by 
N(t) in the space of all matrices. We formulate conditions for quantization of the pumped charge. 



Adiabatic charge pumping has attracted considerable 
attention recently, largely motivated by the experiment 
of Switkes et al.tl A cyclic modulation of gate potentials 
can produce a dc current through an open system, 
contrast to pumps in the Coulomb-blockade regime 
pumping through open systems is governed by quantum 
effects. The pumping frequency exceeds the level spacing 
violating true adiabaticityjj but a compact description 
can still be achieved if the driving fields, change weakly 
on the scattering time scale. BrouwerQ expressed the 
average pumped charge in terms of the loop traversed 
by the scattering matrix S(t). Further research concen- 
trated on the mesescopic fluctuationsj3~l3 the role of dis- 
crete syn^etriesBa and the conditions for charge quan- 
tization,Ell3'E3Eil not necessarily present in open systems. 

Statistics of charge transport are of great importance, 
both for fundamental research and potential applications. 
The counting statistics can be formally expressed^ as the 
determinant of an integral operator involving S(t). While 
compact cpoipssions were obtained for particular pump- 
ing cycleSjtHnU a general case requires further analysis. 

Let us begin by sketching our main results. Our con- 
sideration is valid also in the presence of a voltage bias, 
since the latter can be gauged away at the expense of a 
phase of the scattering matrix. We will show that the 
statistics of current are invariant under a local symmetry 
group and are determined by the path in the correspond- 
ing coset space. Namely, consider the transformation 



S(t) -> U(t)S(t) , U(t) 



U L {t) 
U R {t) 



(1) 



for a scattering matrix S(t) € U(2M). ft causes ad- 
ditional rescattering between M left outgoing channels 
(Ul) and between M right channels (Ur). We show that 
U(t) shifts the probability distribution P(Q) of pumping 
the charge Q by an integer (we set the elementary charge 
e = 1), the relative winding number of the overall phases 
of Ul and Ur: 



p(Q) 



P(Q-W), W=-^lTr(dUU^a 3 ) 



(2) 



Here 173 is the diagonal matrix diag{l_M; — 1m}- I n par- 
ticular, P{Q) is invariant if the path is trivial, W = 0. 



This result implies that the statistics are determined, 
up to an integer offset, by the path N(t) = S^a^S in the 
(matrix realization of the) coset space U(2M)/U(M) x 
U(M). In the single-channel case N(t) = ner reduces to 
a contour C = {n(t)} on the unit sphere in 3D. WeiLnd 
that at low temperature the average pumped chargetJ is 
given by the area enclosed by C on the sphere (cf . Ref . Q6h , 



(Q) — Aiphe 



(3) 



This area is defined only modulo 1 since a surface with 
the edge C can cover any of two complementary pieces of 
the sphere or even cover the sphere several times. A way 
of fixing this integer is discussed below. When the con- 
tour N(t) [but not necessarily S(t)] is small, the pumped 
charge is quantized. Indeed, in Refs. 
was found under these conditions. 
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FIG. 1. A pumping cycle S(t) defines a contour N(t) in the 
coset space, which determines the statistics. The area ^4 sp hero 
enclosed by N(t) in the coset space and the area ^4 m in of the 
minimal surface in the space of all matrices define the average 
pumped charge and the noise. 

The current noise was expressedO via a non-local in- 
tegral involving N(f). We transform it to a local integral 
over the time disk defined as follows: For the analysis of 
driving at frequency to we can replace the time axis by 
a unit circle Ct'. w — e luJt . There is a unique harmonic 
(An = 0) extension of the mapping n(w) from Ct into the 
disk D t . The noise is given by the Dirichlet functional, 



«Q 2 » 



1 

8^ 



(4) 



where w = x\ + 1x2- For a given contour C the de- 
tails of its traversal in time depend on the pulse shape 



1 



of the driving fields. Unlike the average charge (||), the 
noise value is sensitive to the pulse shape, i.e., to the 
time parametrization. Optimal pulse .shapes minimizing 
the noise were found in several cases £3o Here we solve 
the problem of noise optimization for an arbitrary cycle. 
Specifically, we show that the pumping cycle is optimal 
when the mapping uLw) is conforms!, (d w n) 2 = 0. The 
minimal noise valueEj is the area of the minimal surface 
(soap film) spanned by C, 



((Q ))min — . ^.min • 
47T 



(5) 



Similar results hold in the many-channel case. 
Invariance. Physically, multiplication of the scatter- 
ing matrix (Q) by U(t) just redistributes the scattered 
particles, without affecting the correlations at the scat- 
tering center. The outgoing states acquire an extra time- 
dependent phase which changes the time these particles 
need to reach the reservoirs. As a result, the extra charge 
Wi = tfTr(dU}Ui)/4Tri is transferred to the lead i = L 
or R, and we get (g). For periodic U(t) these numbers 
are integer. (Note that no net charge accumulation near 
the scatterer implies Wl — —Wr.) 

More formally, the transformation rule (Q) for (Q) fol- 



lows from the Brouwer formula, 
mulants one can use the resultt. 
function x(A) = ^ Q P(Q)e iXQ : 



For the higher cu- 
for the generating 



x(A) = det l + n F (t',t)(SlM)S x (t)-l) 



(6) 



Here S\(t) = e- lXa3 / 4 S(t)e' Xa ^ 4 and n F (t',t) = 
i/[2ir(t' — t + i0)] is the Fourier transform of the Fermi 
distribution. 



In Ref. 



14 



by separating phases and ampli- 
tudes of S(t), this result was presented in a form, which 
immediately implies the transformation rule (|^) . Indeed, 
the determinant in Eq. (5) of Ref. |l4| is invariant under 
(Q) and the quantity N in the prefactor is shifted by W. 



To express x(A) vm N 

e iXa- 3 /4:^ e -iXcr 3 /2g e iXeT a /4: _ 

Using the identity e ~* AN / 2 = 



we notice that S^^S. 

,iXo 3 /4 

A • • X 



e" 
cos ■ 



iXS f <T 3 S/2 e iX<7 3 /4: 



i sin we get 



X (A) = dot [1 - \n F (t\t)(e lX ^ - l)a a [N(t) - a 3 ]] . (7) 
At T = 0, multiplying by 1 + n F (e~ iA<T3 / 2 - 1), we obtain 



X(A) = dot [l + n F (i',i)(e- lAN( * )/2 - 1) 



(8) 



Note that the result (||) is explicitly invariant under 
global rotations— N(i) — > U^N(i)U [corresponding to 
transformationsO S(t) -> S{t)V}. 

Eqs. (||)~(§) involve N, but not S, and hence can de- 
fine P(Q) only up to an integer offset. Indeed, the infi- 
nite product of the eigenvalues of these integral operators 
can be regularized in many ways. £3 Notice that for the 
operator (g), due to strong degeneracy, one can choose 
eigenstates that span a narrow frequency range, of or- 
der L). For those far above the Fermi level (tip — 0) the 



eigenvalues are 1. Deep in the Fermi sea (rip — 1) the 
eigenvalues appear in pairs e ±jA / 2 with the product 1. 
Although regularization procedures can pair them in dif- 
ferent ways, x(A) can be changed only by an even power 
of e ±lA / 2 , which gives an integer shift of Q. 

Pumped charge and the area. The expression for 
(Q) = <9\x(A = 0) contains a singularity whose regular- 
ization requires the knowledge of the full S(t) and gives 
an integral over the periodtj, 



(Q) 



(9) 



The loop C — {£(£)} can be contracted to a point, 
uniquely up to continuous deformations. In the process, 
the loop spans a surface D. (For a two-parametric pumpEj 
there is a natural choice of D corresponding to the inte- 
rior of the contour in the parameter plane.) Using Stokes' 
theorem, we rewrite (^) as a surface integral, which fur- 
ther reduces to the "area" of the corresponding surface 
V in the coset space, 



(Q) 



Tr{a 3 dSAdS^) 



v 



Tr(NdN ArfN) 
16ni 



(10) 



Note that the integrand is the curvature of the fiber bun- 
dle S — » N = S^a^S. In the single-channel case N = ncr, 
and we obtain (||): 



(Q) = — [ eiikUidn, A dn k . 
8n J v 



(11) 



One can try to define the "integral part" Q lnt of (Q) 
as follows: Let us parametrize scattering matrices as 
S = US°[N] with a matrix S°, defined for any N, and a 
matrix U as in Eq. ([j]), and assign to each cycle S(t) the 
winding number of the corresponding U (t) . This attempt 
fails, since there is no continuous global map 5°[N]. In 
fact, any two loops in U(2M) can be deformed into each 
other, i.e., any <5i nt is discontinuous under certain con- 
tour deformations. However, continuous maps S°[N] do 
exist for contractible regions, and one can introduce Q- m t 
for contours S(t), for which N(t) does not leave such a 
region. Examples are regions of matrices S without per- 
fectly transmitting (or reflecting) channels. In particu- 
lar, the integer N introduced in Ref. |l4| changes abruptly 
when the pumping cycle contains a scattering matrix S(t) 
with a perfectly transmitting channel. 

At this point we can formulate sufficient conditions for 
the quantization of the pumped charge: The fractional 
part of (Q) vanishes for small contours N(t) in the coset 
space, which remain close to their initial point over the 
pumping cycle. An example: the minimal and maximal 
conductance, 5 = and g = M , is achieved at the points 
N = ±03. Hence keeping g close to one of these values 
throughout the cycle guarantees the quantization. We 
evaluate the integral (|10| ) for such cycles and estimate 
the accuracy of quantization as 

SQ < g for g « , SQ < M - g for g « M . (12) 
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For a single channel the fractional part is the area (g) 
within the small contour n(f). 

For a single channel, M = 1, the scattering matrix can 
be parametrized by the conductance g and three phases: 



,03 



The components of the unit vector n then are 
n z = l-2g; n x + in y = -2iy/g(l - g) e i(Q-/3) 



(13) 



(14) 



Using these expressions we can explain the charge quanti- 
zation found, for instance, in Refs. QJn]: For the pumping 
cycles studied the system encircled the resonance point 
g = 1 in the parameter plane at a sufficient distance from 
it so that g ~ throughout the cycle. The corresponding 
loop n(t) encircled the north pole (g = 0; Fig. g). Since 
the interior of the loop in the parameter plane contained 
the resonance point, the surface V in (|ll]) covered the 
lower part of the sphere, i.e., almost the whole sphere, 
(Q)wl. 

Noise optimization. For the noise, given by the A 2 
term i n the Taylor series of ln^(A), we obtain a double 
integrafiaii-J over the unit circles w = e luJt : w' = e lu)t : 



«Q 2 » = 



16tt 2 



dw dw' 
(w — w') 



Tr[l- N(t)N(*')]- (15) 



Let N(t) = J2k>o ^fc exp(ifcwt) + h.c. Then the (unique) 
harmonic extension of the mapping w — ► N(io) from the 
circle into the disk \w\ < 1 is given by 

N(w) = N+(w) + N-{w) = ^2 k>Q N fc w fc + h.c. (16) 



Expanding the integrand (|15|) near t' — t one sees that 
we are justified in replacing the integration over w by 



0. (17) 



w\ — l-\-e 



Eqs. @, (0) and the invariancc under time reversal 
w,w' — > w,w allow us to use complex analysis to do the 
integration over w' in (|l5|): 



<(Q 2 »=Tr 



J 167ri J 



Finally, Stokes' theorem gives the integral over the inte- 
rior Dt = {w — x\ + ix2 : | it; | < 1} of the time circle, 



<(Q 2 » 



l 

167T 



d 2 xTr[(a,N) 2 



(18) 



which reduces t o (|4| ) in the single-channel case. Note that 
the functional (|18|) is well-defined for any surface N(u>). 
The harmonic surface (|l^) provides the minimal value to 
( |l8| ) among the maps N(iu) with the fixed value at the 
boundary, N(t). 



Now we turn to optimization of pumping: The 
cyclic evolution of S(t) is achieved by periodic changes 
in external parameters that control the scattering or 
the bias. For a fixed trajectory in the parame- 
ter space the rate of motion can be varied. These 
changes do not affect the average pumped charge 
( pr| ) but do influence the noise (|l8|). Notice thatE. 
iTr[(9iN) 2 + (9 2 N) 2 ] > VTr[(diN) 2 ] • Tr[(a 2 N) 2 ] > 
area spanned by d{N and 9 2 N in the matrix space. 
Hence the noise value always exceeds the area [defined 
by the scalar product {A, B) = Tr(A+S)/2]: 



((Q 2 » > i-4[NH 



(19) 



The equality is achieved only if <9iN are orthogo- 
nal, Tr(9iN9 2 N) = 0, and have the same length, 
Tr[(diN) 2 ] = Tr[(<9 2 N) 2 ]. These two conditions together 
can be written as Tr[((9,„N) 2 ] = and characterize con- 
formal mappings N(w). 

Since any regular surface can be parametrized con- 
formally, the minimal values of both sides of Eq. jl9| ) 
coincide. The minimal noise value is thus given by 
the minimal area (|^) of a surface spanned by the loop 
N(t). The optimal pumping with this minimal noise 
value for a given loop is achieved when the correspond- 
ing harmonic mapping N(iu) is conformal. Notice that 
the SL 2 (R) time reparametrization symmctrylij, w — > 
(w + a)/(l + dw), preserves the classes of harmonic and 
conformal maps. For N cycles N^r(w) = Ni(w N ) and 
Ni(J| i=1 (?« + a,)/(l + a t w)) give the same statistics. 

Applications. Our findings give a new perspective on 
the analysis of pumping cycles discussed in the literature. 
Consider the cycles 



S a (t) 



_ e ~i<j>(t)a 3 /2 
= p i<t>(t)<?z/2 



S(0} 
S{0)e 



i<f>(t)a 3 /2 



> A(f> = 2tt7V, ( 20 ) 
(21) 



Thp-|fi*st of them, studied extensively by Levitov et 
al.Ji^'tij describes conductors under the voltage bias 
— h<p(t) j e, as one can see by applying a gauge transforma- 
tion. The cycle S a was discussed in Ref. [14[ These cycles 
differ only by a transformation (Q) with U(t) = e^" 3 , 
hence the statistics coincide up to a shift by the wind- 
ing number: Pp(Q — N) = P a (Q). Indeed, the same 
pulse shape 4>(t) = cut [and others, generated by SL 2 (K.)] 
was found optimal for both cycles. The statistics for 
this optimal cycle S a (t) are related to the well-known 
binomial, distribution for a conductor under a constant 
positivM bias by P° pt {N - Q) = P° pt (Q), in agree- 
ment with Ref. |lj. For a single channel the vector n(t) 
follows a line of constant latitude for both S a and Sp: 
(3(t) - (3(0) = cp(t) or a(0) - a(t) = cp(t). The pumped 
charge (Q) is given by the area above this line, g, for 
(^p|), and below this line, 1 — g, for (pi]). The minimal 
noise value @ is the area g(l — g) of the sphere's cross- 
section. The optimal pumping corresponds to the trivial 
homothety of the time disk onto this cross-section, this 
map being obviously harmonic and conformal. 
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FIG. 2. Several pumping cycles n(i): sM and sM lead 
to the same-counting statistics, which at,*= 1/2 coincide with 
that for SjH. The sketched contour SoQ'tJ describes pumping 
of a unit charge per cycle. 

The rotational invariance of xW implies that the 
counting statistics for any-circular path n(i) is the same 
as for biased conductors.tL 2 ] In particular, the optimal 
pumping uniformly traverses the circle and gives rise to a 
binomial distribution. As an example consider the cycle 



s g (t) 



cos i](t) smr](t) 
smr](t) —cosrj(t) 



A?7 = 27r, (22) 



during which the conductance g = sin 2 rj(t) oscillates. We 
find that (Q) = and n(t) traverses twice the meridian 
in Fig. |^. Thus P(Q) coincides with the distribution 
for the equator, Sp at g — 1/2, and the pulse r){t) = 
uit from Ref. [l4| is optimal. For this pulse P(Q) is the 
shifted binomial distribution for N = 2 cycles, P° pt (Q) = 

P 7\Q + 1)1^ 2 = i(q+i)» in agreement with Ref. ||. 

Our geometric approach allows us to obtain relations 
between current and noise for broad classes of pump- 
ing cycles. For small loops n(t) the minimal surface 
lies within the sphere and the analysis simplifies. If 
the loop S(t) is also small, the system is in the weak 
pumping regime with (Q),((Q 2 )) <C 1. In this case 
we have for a general (possibly self-intersecting) loop 
(Q) = (A+ - A^)/4tt < ((Q 2 )) min = (A+ + 4_)/47r, 
where A± are the contributions to the enclosed area 
with positive (resp. negative) orientations. The weak^ 
pumping regime was studied very recently by LevitovtZl 
who found that the transport is described by two uncorre- 
cted Poisson processes (transporting charge to the right 
and to the left), reducing in some cases to a single pro- 
cess. Our inequality ((Q 2 )) > (Q) is in agreement with 
Levitov's findings, with A±/4ir being the rates of the two 
Poisson processes for an optimal cycle. The equality, the 
criterion for the reduction to a single Poisson process in 
the weak-pumping regime, is thus reached only for op- 
timally traversed contours enclosing the area of a con- 
stant orientation (in particular, for non-self-intersecting 
loops). The simplest example of such a cycle considered 
in Ref. [ij] corresponds in our terms to n(t) traversing 
uniformly a small circle. Generally, for weak harmonic 
driving n(i) encircles an ellipse. The optimal pulse shape, 
given by the conformal iHiap of D t onto this ellipse, in- 
volves elliptic integrals! 2 -!! Further, for a general small 
polygon the optimal pumping is given by the Schwarz- 
Christoffcl formula, describing a map of D t onto its into-, 
rior (also reducing to elliptic integrals for a rectangle) .Ell 



The results concerning the weak-pumping regime can 
be generalized to the many-channel case. Using local 
complex coordinates in the coset space, we found that 
the ratio ((Q 2 ))/|(<9)| > 1, reaching the minimal value 
unity (corresponding to a single Poisson processO) only 
for optimal cycles with a complex analytic (or antiana- 
lytic) minimal surface N(u>). 

For the strong-pumping regime our description also 
gives new results. For instance, in the interesting case 
of a single channel and a contour n(t) without self- 
intersection, we find that ((Q 2 )) m in < distance from (Q) 
to closest integer. 

In conclusion, we have shown that the counting statis- 
tics of a sample subject to a periodic pumping (and pos- 
sibly to a voltage bias) are determined by a path in the 
coset space and given by Eq. (||). The average pumped 
charge and its minimal variance for an arbitrary pumping 
cycle are given by the area (10 ) encircled by this path in 
the coset space and the area (5j) of the minimal surface 
spanned by this path, respectively. The optimal pump- 
ing strategy can be found as a harmonic conformal map 
of the time disk onto this minimal surface. Our results 
represent a unifying framework for analysis of transport 
statistics in various realizations of pumping. 
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